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Abstract 

We consider the gravitational collapse of a spherically symmetric ball 
of dust in the general relativistic weak gravity regime. The velocity of the 
matter as viewed by external observers is compared to the velocity gauged 
by local observers. While the comparison in the case of very strong gravity 
is seen to follow the pattern familiar from studies of test particles falling 
towards a concentrated mass, the case of weak gravity is very different. 
The velocity of the dust that is witnessed by external observers is derived 
for the critically open case and is seen to differ markedly from the expecta- 
tions based upon Newtonian gravity theory. Viewed as an idealized model 
for a cluster of galaxies, we find that with the general relativistic veloc- 
ity expression, the higher-than-expected constituent velocities observed 
can be readily correlated with the solely baryonic measure of the mass, 
obviating the need to introduce extraneous dark matter. Hitherto un- 
explained and subject-to-reinterpretation astrophysical phenomena could 
also be considered within this context. It is suggested that an attempt 
be made to formulate an experimental design at smaller scales simulating 
or realizing a collapse with the aim of implementing a new test of general 
relativity. 

Recently pQ, [2J, we demonstrated that general relativity could account for 
the flat galactic rotation curves, the observation of the essentially constant veloc- 
ities of the stars in the galaxies out to their extremities, without the requirement 
for the conventionally demanded vast reservoirs of exotic dark matter. We 
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1 Various critics have claimed that our results stemmed from a singular surface layer of 
mass but we have shown that the singularity actually represents the benign incorporation of 
a discontinuity in density gradient. As well, a different approach [3] largely supported our 
central thesis. 



have seen that general relativistic nonlinearities can play an important role even 
in the context of weak gravity. 

To fortify our claim, and to attack frontally the issue of higher-than-expected 
velocities on the basis of Newtonian gravity and the observed matter in clusters 
of galaxies, we now consider a different problem, the gravitational collapse El of 
a spherically symmetric sphere of dust in the phase where a ball structure is 
evident with external vacuum. While the interest in the collapse problem in 
the past has focused upon strong gravity leading to singularity formation (see, 
e.g. [S]jl, we consider here the weak gravity regime, long before any singularity 
could be formed This could be viewed as a special case model, albeit highly 
idealized, of a cluster of galaxies in evolution where in the normally unsymmct- 
ric case, the component velocities have been observed to have unusually high 
velocities according to Newtonian gravity. Reportedly, this phenomenon was 
the historic origin of the dark matter hypothesis, advanced in the 1930's by 
Zwicky, in an effort to explain the high velocity observations within the con- 
text of Newtonian gravity theory. In this paper, we show that in this idealized 
model, such velocities can be accounted for in principle using general relativity 
in the absence of extraneous assists from dark matter. 

We first consider the treatment in [3] (henceforth referred to as "LL" ) of the 
familiar Schwarzschild solution, the spherically symmetric vacuum gravitational 
field for a spherically symmetric mass m. A spherically symmetric metric can 



In the case of vacuum, the metric functions are readily found by solving the Ein- 
stein field equations and reflecting the intrinsically static nature of the spher- 
ically symmetric vacuum solution, can be expressed in the time-independent 
form 



2 The expression "gravitational collapse" has come to mean the total crunching of matter 
under gravity but we are using it here in the general sense of material flowing inward under 
gravity. 

3 The essential point in that paper was the implementation of pressure during the late stages 
of gravitational collapse. The maintenance of the dust equation of state is totally unrealistic 
physically when the concentration of material mounts dramatically in the approach to zero 
volume. 

4 Part of the motivation for analyzing the present problem stems from the fact that singu- 
larities cannot be raised here as an issue. 

5 When we speak of "high" velocities in this context, we mean high compared to the 
expectation of Newtonian gravity but still very much smaller than the speed of light. Also, 
when we speak of "general relativistic velocity", velocity governed by Einstein's theory of 
gravity, such velocity can take on the whole range of values, but in this paper, the emphasis is 
on the cases where this velocity is much less than the velocity of light. This expression should 
not be confused with the familiar expression "relativistic velocity" of special relativity where 
the reference is to velocity approaching the speed of light. 

6 We choose units where c = G = 1. 
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This is generally referred to as the Schwarzschild metric. LL effect the trans- 
formation to comoving Q synchronous [4] coordinates R,t as d 



t = t+ v ' dr 

1 - 2m 



with f(r) chosen as 



m 

\ (3) 

" = l+ 1 m(i-^) dr 



yielding the simple relationship between the coordinates 

3 -,2/3 



l(R-r)) (2m) 1 /3. (5 ) 

The singularity issues encountered in (TT]) , ([2]) at r = 2m are thus alleviated. The 
metric in these new comoving (_R, r) coordinates is then expressed as 

ds 2 =dr 2 — (hR-r))^ (2m f 3 (d6 2 + sin 2 ed^ 2 ) (6) 

which is now seen to be time-dependent. The strong gravity regime evolves as 
the proper time r approaches R and the singularity emerges at R = r. 

Our focus will be on the weak gravity regime where R >> t for all R. This 
translates to r » 2m for all r in the (r,t) frame. The time coordinate t mea- 
sures time read by the asymptotic observer. The standard general relativistic 
treatment concentrates upon the regime of strong gravity where the difference 
in perception of the proper velocity of a freely falling test particle as measured 
by the local observer in comparison to the measurement of the velocity by the 
asymptotic observer becomes particularly significant (see [1]). 

In generality the proper radial velocity of a freely falling test particle cannot 
be evaluated in the (R, r) coordinates. This is because R is constant for any 
given particle in this frame and hence the radial velocity is always zero in this 
comoving frame. For the required ingredients, LL use the solution of the radial 
geodesic equation for a freely falling test particle in the usual Schwarzschild 
coordinates (r, t) as employed in fl]),© which are suitable for this purpose. 
The geodesic solution for dr/dt and the metric coefficients goo and 511 of |T]) § 
are used to evaluate the proper radial velocity 



-911 dr /yx 

3oo dt 



7 In these coordinates, an observer being at R = constant is in free-fall. 

8 8 and <f> are left unaltered. 

9 (xW 2 ,^ 3 ) = (t,r,0,<t>) 
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This equals y/2m/r in magnitude for particles released from rest at infinity and 
is seen to approach 1, the speed of light, as r approaches 2m. (The rest release 
point tq in LL Eq.(102.7) is taken to be infinite here.) However, for asymptotic 
observers who reckon radial distance and time increments as dr and dt, the 
measured velocity is 

dr ( 2mA [2~^ 

dt=-{ 1 --)^- (8) 

which approaches zero in the very strong gravity regime as r approaches 2m, in 
stark contrast to the proper radial velocity. However, for weak gravity which is 
our focus in this paper, the (1 — 2m/ V) factor in {8| is approximately 1 and the 
local proper and asymptotic measures of velocity are approximately equal in the 
value — y / 2m/r. This justifies the neglect of general relativity in the context of 
weak gravity for the case of motion in spherically symmetric vacuum. 

Interesting new developments ensue when we turn to dust collapse. As 
with the vacuum case, LL choose comoving coordinates for dust collapse and 
structurally as in (J6]), express the metric as 

ds 2 = dr 2 - e^ T ' R) dR 2 - r 2 (r, R){d0 2 + sin 2 0d(p 2 ) (9) 

The four non-vanishing Einstein field equations are extremely non-linear 

-e- A (/) 2 + 2rf + f 2 + l = 0, (10) 

e~ A , „ , , s r\ ■■ X 2 2f , . 

(2r"-r'X') + — + A + — + — = 11 

r r 2 r 

- (2rr" + {r') 2 - rr'X') + ^ + r 2 + lj = 8np (12) 

2(f)' - \r' = (13) 

It is remarkable that for such a complicated non-linear set of equations, at 
least a part of the solution should be as simple as it is: ["1 



f 2 = E(R) + ffl^ (15) 
r 



9F^ ■ 
4 



(R-t) 2 / 3 (16) 



where in (fTB]) . we have chosen the part of the three-stage solution 0] for the 
case where E(R) in (|15p is taken to be zero, the particles being released from 



10 A dot denotes the partial derivative with respect to r and a prime denotes the partial 
derivative with respect to R. 
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E(R) and F(R) are functions of integration. 
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rest at infinity in the infinitely distant past. 00 For all three cases, positive, 
negative or zero E, the density p is incorporated into the solution as 

^P=^2- (17) 

From (fTT|) , a simple integration (see [1]) shows that the mass M(R) within 
the radial coordinate R is 

M(R) = F(R)/2 (18) 

and thus the entire mass M is given by M(Rq) where i?o is the outer comoving 
radial coordinate of the dust ball. 

Our focus here is upon the radial dust velocity measured by distant "rest" 
(i.e. nora-comoving) observers. As with the vacuum case, we must choose new 
coordinates for the evaluation of the radial velocity because the R coordinate is 
constant for any given dust particle. As discussed above, for the dust case, we 
continue to use the approach taken by LL for vacuum and evaluate this radial 
velocity dr/dt in the Schwarzschild-like (r, t) coordinate frame. However, unlike 
the case of vacuum, it is unnecessary to solve the geodesic equations at this point 
because the motion of the dust medium has already been solved in the comoving 
frame. It is this motion that is of concern to us£j What is required is to re- 
express the solution in Schwarzschild-like (r, t) coordinates. For consistency 
with the solution form of (|16[) and to maintain maximum available generality, 
we choose the general form of transformation with arbitrary functions p(r, t) 
and q(r, t), 

s/FR = p(r,t),VFT = q(r,t) (19) 

with the constraint 

p(r,t)-q(r,t) = (2/3y/ 2 . (20) 

12 There are two other cases for non-zero E, the bound and the unbound cases, which are 
familiar in concept from classical mechanics. These solutions are expressed in parametrized 
form in [4]. Here, we focus upon the simplest critically open case. 

13 Note that in general, one can choose an arbitrary function tq(R) in place of R in the 
(R — t) 2 / 3 factor in l|16J l, We have chosen to(R) = R to mesh smoothly with the chosen 
vacuum Schwarzschild solution in comoving coordinates j5)which holds in the vacuum region 
surrounding the ball of collapsing dust. 

14 Note that there is no problem in matching the interior dust solution with that of the 
exterior vacuum Schwarzschild solution and no reasonable argument can be raised as to the 
presence of a surface layer. A well-known example of interior-to-exterior matching with spher- 
ical symmetry is that of Schwarzschild's static constant density ball with pressure matched to 
the exterior vacuum solution. 

15 Actually, the motion is also geodesic here since there is no pressure. 

16 Within the dust ball, the coordinates used cannot be expressly Schwarzschild coordinates 
as used with explicit time-independence in (|2}because the field within the dust is intrinsically 
dynamic. However, the essential metric structure is the same. As well, by making r 2 the 
coefficient of the angular part of the metric, the circumference of a ring of particles at r 
assumes the familiar flat-space value 2nr for both the proper measure and for the measure as 
judged by distant observers. 
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From |gn|), we see that0 

p'(r,t) - q'(r,t) = r^ 2 ,p(r,t) = q(r,t). (21) 

We take differentials of ([11]), and with (|2"0 )) .([2"T j) . solve for dR and dr. These 
differentials are substituted into ^ to derive the normal form of the metric in 
Schwarzschild-like coordinates (r, t) with terms of the form goodt 2 and g rr dr 2 , 
as well as an undesired cross-term of the form 2go r drdt. This cross-term must 
vanish to mesh with the exterior Schwarzschild metric at the vacuum interface 
and maintain the useful Schwarzschild-like non-rotating form within the ball. 
This metric form includes a yet-to-be-determined p'(r, t) which we set by making 
ffor = yielding 



P 



where 



(a + /3)(l-/? 2 ) 

_ rF' _ rM'(R) 
~ 3F ~ 3M{R) 
j= jF_j2M{R) 



(22) 



(23) 



Also required in the calculation for p' is e x which, from (jT4j) , is equal to (r') 2 
for E — 0. In turn, this requires r' which is computed from (I1G|) yielding 

r' = a + /3. (24) 

Since the R coordinate is comoving with the matter, we express the condition 
for the radial motion of the particles by taking differentials of the first of (|19p 
and setting dR = 0: 

p'(r,t)dr +p(r,t)dt = (25) 

from which we find the form of the radial velocity of the particles as witnessed 
by external observers 

dr/dt = -p(r, t) /p'(r, t). (26) 
To solve for p(r,t), we first apply d/dt to (fTT|) : 

St r J( Q + / j)J(M|E> 



17 Here and in any subsequent appearances of p(r, t) and q(r,t), a dot on these functions 
denotes the partial derivative with respect to t and a prime on these functions denotes the 
partial derivative with respect to r. 

18 For motion in this spherically symmetric study, d0=d<f>=0 as well. 
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The derivation of (J27J) made use of 



d(a + (3) 
dt 



F' r F" 
2y/Fr o~F~ 



p>2 1 
T 2 



OR 
~dt 



(where (|23|) has been used) and the elimination of 9/l 



dt usm S 



2s/F 



OR 
dt 



(28) 



(29) 



which follows from the partial differentiation with respect to t of the first of (fl9|) . 
Finally, using (|22| and( |27|) in conjunction with (|26|) (and with a cancellation of 
the factor ( 3R ^f a + y/ffi)), we find0 



dr _ (a + /3){l-P 2 ) 
dt ~ 87rrV 

This is the key equation. The complexity of this velocity expression as 
computed by observers external to the distribution of matter is in very sharp 

contrast to the simplicity of the proper velocity form (3 — \J~^r as witnessed 
by local observers. However, it is the former that is relevant for astronomical 
observers. For local observers, it is the mass, F(R)/2 = M(R), at radii within 
the point of interest that determines the velocity, the situation as in Newtonian 
gravity. However, we see in (|30[) that within the context of general relativity, the 
external observers ascribe velocity based on additional factors, the reciprocal of 
the local density squared and its time rate of change (also expressible as the 
time rate of change of reciprocal density), the gradient of the mass within the 
radius in question, M'(R) as well as well as its gradient, M"(R). 

We also see that in the limit of very strong gravity, with (3 approaching 1, 
the situation is the same as we witnessed in vacuum: the local observers see the 
velocity approach 1 whereas the external observers see the velocity approach 0. 

However for weak gravity with /3 « 1, the vacuum and dust comparison is 
very different. While the local and asymptotic velocity measures for observers 
plotting freely falling test particles in vacuum in the field of a concentrated 
mass are approximately the same, namely (3, the corresponding velocities for 
local and asymptotic measure in the case of dust are very different in general: 
the velocity is simply (3 for the local measure whereas the asymptotic measure 
is given by the rich expression (|30|) with 1 — (3 2 approximated by 1. Indeed, 
given the complexity of the form of dr/dt in (|3"0"|) , it would be a very special 
occurrence for dr/dt to have the value (3. Thus, when astronomers witnessed 
with consternation velocities greater than (3 in galactic clusters, they should 
have been more reasonably surprised had they witnessed (3 velocities. 

19 Note that a lengthy calculation with the metric components in the (r, t) frame in con- 
junction with J7J yields the proper radial velocity, the velocity measured by local observers, 

with value \ — . This coincides with r = However for observations by external observers, 

y r or J 

the measured velocity is dr/dt as given by H30H 4 . 



a F" 
F \ (F') 2 



1 

2F 



2E 

dt 



(30) 
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As an application of (f3"0")l . we first consider the astrophysical realm. Most of 
the gravity in the universe is weak gravity where Newtonian theory has been 
deemed to be perfectly adequate. Thus, when galaxies in clusters, with gravity 
found to be weak, were seen to have velocities exceeding /3, dark matter was 
introduced as the necessary mass booster to align the observations with enlarged 
P. Newtonian theory formed the basis for the calculations. However, we have 
seen that general relativity, which is almost universally accepted as the preferred 
theory of gravity, actually predicts velocities that have elements beyond (3 even 
when the gravity is weak. A key point is that the nonlinearities of general 
relativity play an important role in this problem, leading to expressly non- 
Newtonian behaviour, even though the gravity is weak. 

For the Coma Cluster of galaxies, the ratio3 2M(R )/r is of order 10~ 4 
assuming the existence of dark matter and of the order 10~ 5 by not assum- 
ing any dark matter. The gravity is indeed very weak in this source for the 
kind of applications under consideration. Thus as a test model, we consider an 
idealized Coma Cluster of galaxies, one of spherical symmetry with the veloci- 
ties as reported in [Sj. At a radius of 1 Mpc, the total cluster mass, including 
dark matter, is given as 6.2 x 10 14 M Q , with the 13%-17% portion being normal 
baryonic matter. Within a radius of 3 Mpc, the total mass is reported to be 
1.3 x 10 15 Mq, with the normal luminous matter portion within the wide range 
of20%-40%. 

It is simple to fit these data with an accumulated mass function 

F(R) = kiR k2 , (31) 

(fci, k2 constants) as shown in Figure [TJ With F(0) = from (|3ip . we are 
assured that there is no singularity at the origin 0]. Using (|3"T|) in (|17p . we 
derive the density profile for the distribution. The graph of the densities for the 
two extremes of the uncertainty range and the average is shown in Figure [2J 

The velocity associated with each F(R) is given by (f3T)]) where we can set 
the "boosted" velocity as 

dr/dt = -n/3 (32) 

where j3, as throughout the paper, is composed from the baryonic mass alone 
and n is the "booster" number to bring dr/dt to the observed level of velocity. 
Assuming the baryonic mass is 20%, 30% and 40% of 1.3 x 1O 15 M , we find 
that the boost factors n are 2.23, 1.82 and 1.58, respectively. Applying this 
to (j3"0|) , we can solve for dp/dt, the sole unknown factor. The results are: 
2.13 x 10- 41 kg/m 3 /sec, 2.62 x 10" 41 kg/m 3 /scc and 3.02 x 10" 41 kg/m 3 /sec, 

20 The key equation (1 30 P> is most conveniently expressed in terms of r and _R. While r has a 
direct measurable connection to the source in that Anr^ is the surface area of the ball, there 
is a great deal of arbitrariness attached to the choice of the numerical value of Rq. From 
the transformation equations, we see that different settings of zero value for the clocks will 
change the number attached to Ro- For the purpose of normalization to enable numerical 
calculations, we must choose Ro = ro so that the average density of the ball calculated using 
ro will equal that using Ro. The zero setting of the clock is adjusted to assure this equality. 
With the average densities normalized, the explicit calculation of dr/dt can proceed. 
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Accumulated Mass vs Radius 
6e+14 i 




Mpc 



Figure 1: The upper, middle and lower limits of mass accumulation are described 
by the functions, F = 6.641 x lO" 16 ^' 453 , F = 1.244 x lO" 12 ^ 305 and F = 
2.531 x lO" 7 ^ 066 respectively. 
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respectively. Rates such as 10 _41 kg/m 3 /sec are quite reasonable as over a period 
of one billion years, the density would grow by 10 -25 kg/m 3 , hence roughly 
doubling the value of the present density. 

In this example, we see adequate scope to explain the observed velocities 
within the framework of general relativity without the requirement of any ex- 
traneous dark matter. The new elements of local density, its time rate of change, 
the gradient of the mass interior to the observation point as well as its gradi- 
ent are additional factors that ultimately determine the net observed velocity 
of the matter by external observers. While this is an idealized case of per- 
fect spherical symmetry, it would seem reasonable to expect comparable effects 
for non-spherical accumulations of freely-gravitating collections of bodies as we 
have in clusters of galaxies. Had Zwicky made this calculation 70 years ago, he 
might have come to very different conclusions. 

Clearly there is considerable further analysis ahead. This paper has only 
dealt with the simplest case E(R) — 0. The positive and negative cases for 
E offer greater freedom of expression. Ultimately, the ideal would be to for- 
mulate the equivalent effects of general relativity as applied to chaotic weakly- 
gravitating systems. For this, the general relativistic equivalent of the virial 
theorem is called for. As well, there is the issue of the interpretation of lensing 
as a mechanism for the deduction of mass. The subtleties of general relativistic 
weak gravity that we have found in the present work must now be directed to 
the consideration of lensing. 

It must be stressed that as before [T], [2], we are witnessing here the power 
of the nonlinearities inherent in general relativity in the context of weak gravity 
to effect very significant changes relative to the results expected on the basis 
of Newtonian theory. It suggests that hitherto unexplained astrophysical phe- 
nomena be re-considered on the basis of the application of general relativity to 
weak-field gravity. Indeed the present case also suggests that an attempt be 
made to formulate an experimental design at smaller scales simulating or real- 
izing a collapse with the aim of implementing a new test of general relativity. 
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